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Abstract We develop a quantum replica method for interacting particle systems 
and use it to estimate the location of the glass transition line in Helium 4. Although 

C^ ■ we do not fully succeed in taking into account all quantum effects, we make a 

thorough semiclassical analysis. We confirm previous suggestions that quantum 
fluctuations promote the formation of the glass and give a quantitative estimate of 
' O I this effect at high density. Finally, we discuss the difficulties that are met when one 

C , tries to extend the calculation to the region of low densities and low temperatures, 

where quantum effects are strong and the semiclassical expansion breaks down. 
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1 Introduction 
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^n 

t^^ ■ A series of recent experiments on solid He^, initiated by Kim and Chan [1,2] 

CN I and then performed by several groups [3], raised, among many others, the impor- 

f>^ ■ tant question of whether disorder and fast quenches can induce the formation of 

l^ I a quantum glass phase in Helium 4 and to what extent this new phase is related 

to the super-solid behavior observed experimentally. Following early numerical 
and theoretical analyses, that showed that perfect crystals do not support super- 
solidity [4, 5], and a numerical analysis that proposed the presence of a super- 
glass phase at intermediate densities [6], Ritner and Reppy [7, 8] showed that 
^% I fast quenches produce disordered samples with a change in the moment of inertia 

that corresponds to an extremely high superfluid fraction, of the order of 20%. 
Since then, whether this behavior can be attributed to the presence of a quantum 
(super-)glass phase has been the subject of some debate [9, 10, 11, 12, 13, 14, 15]. 
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However, before addressing the difficult question of whether supersolidity is pos- 
sible, present and possibly enhanced in quantum glasses, it would be desirable to 
have a full understanding of the nature of the glass transition in presence of strong 
quantum fluctuations, a subject which is instead very poorly understood. 

Indeed, several authors attempted to build a theory of the quantum glass tran- 
sition [16, 17, 18, 19, 20, 21], but none has so far been able to obtain quantitative 
results for realistic systems such as Helium 4. One of the main difficulties is that 
new tools are needed to treat glassiness for quantum interacting particles systems 
in a consistent way (the majority of previous works focused on spin systems and 
coarse-grained field theories). Recent mean field studies based on Quantum Mode 
Coupling Theory (QMCT) [22] and on lattice models [23] gave new insights, and 
a coherent qualitative picture of the glass transition in quantum hard spheres has 
now emerged [23, 24]. The result is somehow surprising. Basic intuition suggests 
that increasing the strength of quantum fluctuations should enable atoms to better 
explore the phase space through tunneling, and thereby inhibit the onset of dynam- 
ical arrest that is necessary for glass formation. In contrast, QMCT predicts [22] 
that increasing quantum fluctuations, by augmenting the thermal wavelength of 
the particles (e.g. by reducing their mass), actually favors the glass phase. Indeed, 
one finds that the glass phase emerges at lower density. This result has also been 
obtained in the exact solution of mean field quantum lattice glass models [23], and 
is corroborated by path integral numerical simulations [22]. The analysis of a toy 
model provided the following intuitive explanation [25]: The classical phase space 
of a dense particle system is split into many different basins; in the liquid phase the 
system visits a lot of basins that are not very efficiently packed and therefore have 
a smaller free volume (the volume accessible to each particle to vibrate around 
its equilibrium position). In the glass phase, the system instead visits a few rare 
basins that are very efficiently packed and therefore have a larger free volume. 
Adding quantum fluctuations, particles can lower their kinetic energy by delo- 
calizing, therefore the kinetic energy gain is proportional to free volume: hence, 
glassy states have lower kinetic energy and are favored by quantum fluctuations. 
The net result is that the glass transition line is re-entrant, moving at lower density 
on increasing quantum fluctuations [22]. For strongly quantum systems, however, 
the transition changes nature and becomes a first order phase transition between 
a superfluid phase and the glass [23]. This first order transition is accompanied 
by phase coexistence, which leads to an heterogeneous "superglass", where re- 
gions of low-density superfluid liquid coexist with regions of a high density glass. 
See [23, 24] for a more detailed discussion. 

The results mentioned above have been obtained on simplified models, there- 
fore they only provide qualitative guideline for the physical behavior of Helium 
4. Quantitative analytic calculations of the phase boundaries for systems such as 
Helium 4 are still missing. Unfortunately, Mode-Coupling Theory is of limited 
use since it is known to give poor estimates of the transition point and, at least at 
the present stage, neglects exchange effects [22] . The aim of this paper is to adapt 
the replica theory of glasses [26, 27, 28, 29, 30], which has been shown to give 
good quantitative estimates of the glass transition point for classical liquids, to the 
quantum case; and then apply it to Helium 4. 

The paper is organized as follows. First, in section 2, we use the standard 
prescription of replica theory to derive an expression of the glass free energy in the 



semiclassical region. We show that this is enough to obtain a quantitative estimate 
of the re-entrance effect at high density. In section 3 we speculate on the possible 
phase diagram of Helium 4 in the high density amorphous metastable region. Next, 
in section 4, we discuss the assumptions we made and identify the main problems 
of the theory: we show that there are some important conceptual difficulties that 
prevent us from extending the calculation to the strongly quantum regime and from 
including exchange. While we believe that these difficulties can be in principle 
overcome, we leave this task for future work. 



2 The semiclassical replica method 

We start our discussion by a straighforward application of the classical prescrip- 
tion of the replica method to the quantum case. We consider a (i-dimensional sys- 
tem of quantum Bosonic particles, of mass M, interacting via a potential v(r). As 
usual h is the Planck's constant and h = h/2n. We set the Boltzmann constant 
ks ~ I. If Xi is the position of particle ;', with ;' = 1, • • • ,A'^, the Hamiltonian is 

where Vxj is the gradient with respect to jc,. 

In the following we will make an extensive use of the replica theory of classical 
glasses. We refer to the original works [26, 28] and the review [29] for details. 
The prescription of the theory is to consider a system made of m copies of the 
original system, subjected to an inter-replica attractive coupling that we choose 
to be harmonic [28]. We denote by Xai, a = I, • • • ,m, the copies of particle ;'. The 
Hamiltonian of such replicated system is 

" = -if^I^l + I K^«< -^«y) + ^ I i^ai-XMf , (2) 

^^^' ai a.i<j ^'" i.a<h 

where, compared to [28], e has been rescaled by m for later convenience. Since we 
want to investigate the semiclassical limit, we neglect exchange for the moment, 
and consider the particles as distinguishable. 



2.1 Small cage expansion 

For each molecule ;' (we now omit the index i) we consider the column vector x 
such that x^ = {xi,- ■ ■ ,Xm)- We perform the change of variables ya = v^x with 
v[ = -^(I,- •• , I), v^vi = and v^v^ = 5ab- The Jacobian of the transformation 
is 1 because the vectors Va are orthonormal. Thus the kinetic term is invariant 
under the transformation, Y.a '^^ ~ ^a ^y^ • Defining the matrix of all unit entries, 
lab = 1^ denoting by 5 the identity matrix, and using the completeness relation 
5 = LaVavl, we can write La<h{xa - Xhf = x'^{m8 - l)x = tahx'^Vavl{m5 - 
i)vbvlx = Y^ahyavl {m5 — i)vbyh ■ It is easy to show that {m5 — I)vi = and 



that {m5 — l)va = niVa for a> I. Then, Y,a<h{^a — Xj,)^ = YJa^i'^yl ■ Thus the 
Hamikonian in the y variables reads 

^ = -|^IVL+ I v(^.-^./) + | I fai, (3) 

ai a,i<j i,a>l 



where for each molecule ;' 

yi 






(4) 



with Z = y\/ \/m = (La^a)/'" the center of mass of the molecule. 

Now, following [28], we perform a large e expansion. For large e we assume 
that the ya, a > 1 are small due to the harmonic term, and we expand the potential 
v{xai —Xaj). For a given pair of indices i,j (that we now omit), we define Axa = 
Xai~Xaj, and it follows from Eq. (4) that Zijc^ =ZiZ + J^^>](vft)aZi}'^. We denote by 
jU, V = 1, • • • ,d the spatial indices, and 5^v = dv/dx^ the derivative with respect 
to coordinate x^. Then we can easily expand the potential as follows: 



Y,v{Axa)=mv{AX) + Y,d^iV{AX)Y,ivb)aAy'i 

a a h>\ 

+ Wd,,v{AX)l^{v,)Ay>l^l^{v,)Ay: ■ 



(5) 



^ a b>\ Ol 



Recalling that for b>\, Y.aivb)a = \/mv\vb = and that Y.a{vb)a{vc)a = vlvc = 
5},c, we get Y.aV{Axa) = mv{AX) + \df,yv{AX)Y.b~^^Ay'^Ayl . Collecting the 
terms together and changing variables to Z, = yuj \fm, we finally get 






2M y-a. 2^^"'' 



Ho{X) ^harm(3'^;) 

■^ i<j a>\ 



(6) 



At this order, the Hamiltonian is the sum of a term Ho {X) that describes the system 
of the centers of mass, with a renormalized mass mM and renormalized interaction 
mv{r), plus a sum of independent terms H\^^^^{y), that describe harmonic oscilla- 
tors of spring constant e. The centers of mass are coupled to vibrations by the term 
Vint; since the vibrations are small for large e, this coupling is small and we can 
treat it as a perturbation. 



2.2 Computation of the partition function 



We now want to compute the partition function of the system at temperature T, 
using the Hamiltonian in (6). Again, we stress that in the following we consider 
that the particles are distinguishable. Otherwise, even in the absence of the cou- 
pling Vint, one should consider the y variables as "internal" degrees of freedom 
(like spins) of the centers of mass, and one should therefore impose the permu- 
tation symmetry on the global wavefunction, that would correlate the X and y 
parts of the wavefunctions. Instead, for distinguishable particles, there is no spe- 
cial symmetry requirement and we consider that in absence of Vint the space of 
wavefunctions is a product of independent wavefunctions for the X and y parts. 
Under these assumptions, we get, developing at first order in Vint: 



Z,„(e)=Tre-'^^ = Trxe 



-PHo(X) 



Or. / \\ Nd(m—\) 



where (Vint) is the quantum and thermal average of the perturbation over the X 
and y. The partition function and free energy of each independent oscillator is 



Zharm(e)=Tr,e-^^harm(v 



(y) 



g/3/ifl) _ I 



^harm(e) = -7'logZharm(e) = -T'log 






(7) 



■ t^harm(£) - TSharmi^) ) 



where the frequency (0 = y^e/M; there are d{m— 1) independent oscillators per 
molecule. It will also be useful to define the mean square displacement of the 
oscillator: 



{i 



de 



hO) 
e 



1 



1 



2 e/S»«o - 1 



— n[(0) 
e 



(8) 



with n{(0) = 1/2+ I/(e^'*™ - 1) . 

The average of the interaction term over y is 



(Vmt), = ^l,d^vv{X,-Xj)l^{iy,^yj)^{y,~yj):)y 

i<j a>\ 



{m-\){y^),Y,Av{Xi-Xi), 



(9) 



using that different spatial components and sites i and j are uncorrected, ( (y,- — 

yj)a{yi-y,iYa)y = ^iiviiy^af + {y%f)y = 25f,v{y^)e- Taking the average of Vint 
also on the centers of mass, we get the first order correction to the free energy: 



<P{m,e) = --logZ,„{e) 



-Td{m— l)log 



e/3fiffl/2 



- ^ logTr;,.-/^«<'(^) + (m - 1 )(/),( 1 £ 4 v(Z,. - Z,.) 



(10) 



A^.' 



=c/(m- I)Fhann(e)+^oH + (m- I)(/),^ / drAv{r)go{r) 



where FQ{m) = — T/A/^logTrxe"^^"'^' is the free energy of the effective liquid 
of the centers of mass, and p^go{X -Y) = {£,1^1 S{X -Xi)5{Y - Yj))q is its 
pair correlation function. In the following we will use the notation {Av{r)) = 

^JdrAv{r)go{r). 

Note that the Hamiltonian//o(^) can be rewritten as //o(^) =ot[ — L; J^Tm^x' 
Y,i<j v{Xi — Z,)] , then the fluid of centers of mass can be viewed as the original 
fluid at inverse temperature j5 — m/3 and with a mass M = m^M. Therefore: 

-Fo(m) = -i-logTrxe-'^^of^) =F^^{m^yM) , (11) 

m niN 

where Fiiq(j3,M) is the free energy associated to the original Hamiltonian (1). 

2.3 The replicated free energy 

Recalling that ^^S^hiM = i(/)e, we rewrite (10) as 



de 

^2» 



4>{m,e) = mFiiq{mP ,m M) + d{m - 1) 



p , . , ^il>arm(e) 2(Ziv(r)) 



de d 



(12) 



mFiiq{ml5,m'-M)+d{m-l)Fh.^,^(e + K) 



having defined the effective spring constant K = 2{Av{r))/d. This is correct as 
long as k: ^ e, which is indeed the case since we are working in the limit of large 
e. Now, following [28], we have to extrapolate the free energy obtained in the 
large e expansion to the limit of zero coupling e. This can be done via a Legendre 
transformation with respect to e, following closely the classical derivation of [28]. 
This procedure provides a formally correct treatment that shows that the large e 
expansion is equivalent to a small cage expansion [28]. However, we see that in 
Eq. (12) can set directly e = and we obtain a meaningful result for the repli- 
cated free energy, which indeed coincides with the one obtained via the Legendre 
transform procedure. For the sake of coinciseness, we adopt this procedure here: 
see [28] for further details on the small cage expansion procedure. Finally, we get 
for the repicated free energy: 

<&(m) = mFiiq(m/3,m2M) + ^(m- l)Fharm('c) , (13) 

where the function fharm(£) is defined in (7) and k = 2(ziv(r))iiq/ii, where the 
average is over the liquid with temperature T /m and mass trP'M. 

Given <P{m), the thermodynamic properties of the glass, i.e. the internal en- 
ergy /(j8,m) and the complexity Z(j3,m), are obtained following [26, 28]: 

//3,m=^^^, Z{P,m)=Pm' ^^^ . (14) 

am dm 

Using the relations S\i^{T) = -dTFiiq{T,M) (entropy), Kn^iT) = -MdMFxi^{T,M) 
(kinetic energy), we get for the equilibrium complexity (i.e. the complexity at 
m = 1): 

L,^(T) = LiP, 1) = 5iiq(r) - 2l5Kii^{T)+l5dFi,,,UK) 

= Siiq{T)-dSi,,,UK)-li[2Kii^{T)-dUi,,UK)]. 



In the classical case the last term vanishes due to equipartition, since [/harm = T 
and A'liq = dT /2. This is not a priori true in the quantum regime, but it should be 
true if we assume that the liquid is a superposition of "glassy" metastable states 
corresponding to harmonic vibrations around amorphous positions. We will com- 
ment later on this crucial issue. 



2.4 The classical limit 

Let us first consider the classical limit fi^O. We have n(ft)) -^ (j3/j(jt))~\ and 

iliarm(e) ^ T\og^h(0 = TlogA + - log -^ , (16) 

z zni 



where A = h/\j2nMT is the thermal wavelength. Moreover the liquid of the cen- 
ters of mass is classical and has a partition function 

Z^q = j fl^e-l^^- i-P'nL<A^i-Xj) =A-''^m"''I^Z,{fim) , (17) 

where Zc(j3) = {Niy^ Jd^Xe'l^^xj"^^''^!^ is the configurational partition func- 
tion. Therefore, Fo(ot) = mFiiq(ml5 ,m^M) = t/ TlogA — ^logm+mFc(j3m), where 
FciP) = -{T/N)logZc{P). Using these results, Eqs. (13) and (15) become 

<P{m)=dmTlogA-—-logm + mFc{l5m)+ ^ ' log ^ ^ " 



2 ° '^'" ' 2 ° %Td 

d (Mrrr- 

2 ^ nTd 



L,^iT)=Si,^iT) + dlogA + ^log^-^^^^-d (18) 



where Sc{T) = 5'iiq(r) +dlogA — d/2 and (•)* denotes an average over the clas- 
sical liquid at temperature T* = T /m (with m = 1 in the second line). These are 
exactly the results obtained by Mezard and Parisi in the classical case [28]. 



2.5 First order quantum correction 

We now compute the first order quantum correction to the complexity in the semi- 
classical expansion. This can be done using a formal expansion in powers of 
h [31]. For simplicity we will perform the computation directly for m = 1, starting 
from Eq. (15). The correction to the free energy of the liquid is 

5^iiq = 4^(4v(r))*, (19) 



then it is easy to see that 

a2 n 1 

(20) 



d[Su^{T)-2l5K,,^{T)] = -^ 



^{Av{r)r+dT{Av{r)y 




T(K) 




T(K) 



Fig. 1 The complexity Es;q{T), both in the classical limit, Eq. (18), and with the quantum cor- 
rection, Eq. (23), for two representative densities as a function of temperature. 



To compute the variation of the free energy of the oscillator, we must take into 
account that the spring constant K also has quantum corrections. From (19) we get 



^So{r) = A^Agoir), and using this we get 



f«. = f;*< 



Jdr5gQ{r)Av{r) = ^ J dr^^Ago(r)Av{r) = — (zi2v(r))* . 



(21) 

Denoting Kq = 2{Av{r))*/d and (Oq = y^Ko/M the classical spring constant and 
frequency, the first correction to the free energy of the oscillator is: 

cro.^ . M Aphcoaf d5K a2 , ^ ^_ dA^{A^v(r))* ^^^^ 
8[l5dFk,,^{K)]=d'-i^^^ + -— = -—-{Av{r))* + ^^\^44/^. (22) 



24 2 Kb 24nT 

Collecting these results we obtain 

A2 



4871 {Av{r)Y 



dL{T) = 



4871 



:''^ -*<-"')>■ 



(23) 



Note that using the results above one can check that at this order in h we still have 
2^iiq = '^t/harm('^) SO the last term in (15) is still zero. 



2.6 Results 



We now present and discuss the results of the semiclassical expansion presented 
above. Note that altough Eq. (23) has been derived in a formal expansion in pow- 
ers of h, we will in the following apply it using the real values of all physical con- 
stants. We focus on Helium 4, which we describe by means of a Lennard-Jones 

(LJ) potential v(r) = 4e [(f) ^^ - (f )^] with e = 10.22 K and a = 2.556 A [32]. 

Given the Helium 4 mass, we have h /(2M) = 6.0596 K A . The intensity of 
quantum fluctuations can be conveniently quantified by the ratio of the thermal 

wavelength and the inter-particle distance, A* =Ap^/^ = p^l^JAnh^/{2MT) — 



8.726AK^/^ X pV37'-i/2 Although better potentials have been constructed [32], 
we choose the Lennard-Jones because it is the simplest and it gives reasonable 
quantitative results; in addition, the glass transition of a classical monoatomic 
Lennard-Jones system has been studied in great detail both analytically [28] and 
numerically [33], so we will be able to compare the results obtained in the classical 
case with the first order quantum correction. 

To compute the classical equilibrium complexity and its quantum correction, 
according to Eqs. (18) and (23), we need the entropy of the classical liquid, as well 
as its g{r), as a function of temperature. Note that since m = 1, the averages (•)* 
reduce to standard thermal averages over the classical liquid, that can be computed 
using g{r). These quantities can be obtained using one of the integral equations 
of classical liquid theory [34]. Following [28], we choose the HNC closure [34] 
which has the advantage of providing a direct route to compute the free energy, 
and gives thermodynamic quantities with reasonable accuracy (of the order of 
10%). The HNC equations for the LJ potential have been solved numerically using 
standard iteration schemes. 

In Fig. 1 we report the complexity Zeq(7^) (also called configurational entropy 
in the context of structural glasses). The point where it vanishes marks the Kauz- 
mann transition Tk to the ideal glass. We see that although the quantum correction 
to Zeq(r) can be positive or negative, close to Tk it is always negative, therefore 
increasing the value of 7k with respect to the classical limit. At the lowest densi- 
ties, Tk is found to increase by a factor of two with respect to the classical value. 
It is worth to note, however, that the semiclassical expansion is particularly bad 

o -3 

for p < 0.055 A : in this region, the classical Lennard-Jones system enters the 
gas-liquid phase separation region [35], and the HNC equations cease to converge. 
Instead, for quantum Helium 4 the phase coexistence region is known to be much 

o -3 

smaller, ending at p ^^ 0.022 A [36]. This is due to the large zero-point kinetic 
energy, and it signals that a semiclassical expansion is definitely not reasonable at 
low density, as expected. 

In Fig. 2 (left panel) we report the glass transition temperature 7k (p), obtained 
from the classical and the quantum computation. We see that at all densities, the 
quantum value is larger than the classical one (with the relative difference decreas- 
ing on increasing p). The classical value compares well with the value obtained 
from a numerical simulation of the classical LJ fluid [33]. The quantum value is 
below the melting temperature of Helium 4, as it should. Unfortunately, as already 
discussed, we cannot obtain results for the really interesting region of low den- 
sity and temperature, where quantum effects are most important. For the classical 
LJ potential, Tk increases with density, therefore the transition line in the {p,A) 
plane is re-entrant. This is a trivial classical effect and has nothing to do with the 
re-entrance found in [22, 25, 23]. The important point, for the LJ case, is that 
quantum corrections increase the value of 7k, thereby promoting glass formation 
as stated in [22, 25, 23]. The region A* ^ 1, where the semiclassical computation 
breaks down, is where we expect that the system should undergo a phase transition 
towards a superfluid state (the line A* = 1 is reported in Fig. 2 for comparison). 
Since exchange was neglected from the very beginning, we cannot obviously an- 
alyze this transition. 
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Fig. 2 (Color online) (Left panel) Glass transition lines Ti^{p), defined by £eq(7") = 0, as ob- 
tained from the classical limit Eq. (18) [28] and including the first semiclassical quantum cor- 
rection Eq. (23). Also shown is a numerical determination of the classical glass transition Tg{p) 
from Ref. [33], and the numerically and experimentally determined melting temperature rni(p) 
of Helium 4, obtained from Refs. [36, 37]. The line A* = Ap^^^ = 1 below which quantum 
effects are very strong is plotted as a dashed line. (Right panel) Complete temperature-density 
phase diagram of Helium 4. Equilibrium phases (from Refs. [36, 37]): The red line delimits the 
liquid (L)-gas (Ga) coexistence region. The green line is the A -line that separates the liquid from 
the superfluid (SF). The blue line is the melting transition below which the crystal phase is stable, 
and the disordered phase is metastable. Glass phase: The purple line is the glass (G) transition 
line obtained in the present work within the first order semiclassical expansion. The behavior 
of the amorphous phase, say, within the dashed black line is currently not well understood. The 
two black dots mark the state point investigated numerically in Ref. [6]. 



3 Speculations on the phase diagram of high density metastable Helium 4 



In Fig. 2 (right panel) we report in the density-temperature plane^ all the infor- 
mation now available for the phase diagram of Helium 4. We also draw the glass 
transition line that can be reached if crystallization can be avoided, see e.g. [38] 
for an experimental investigation of supercooled liquid Helium. Ref. [6] reports 
the existence of a metastable disordered superfluid state at T = 0.2 K and densities 

o -3 o -3 

p = 0.0292 A and p = 0.0359 A . These state points are marked by a black dot 
in Fig. 2. However, it is not clear to us whether this state is a liquid or a glass^. 

' The reason why we do not discuss the phase diagram in the more common pressure- 
temperature plane is that the computation of the pressure from the semiclassical expansion 
would introduce additional uncertainty on the position of the glass transition line. Moreover, 
both theoretical and numerical computations are more naturally done at constant density. 

^ In Ref. [6] the system was named a glass because it was observed that a non-homogeneous 
density profile lasted for a certain time during a Path Integral Monte Carlo (PIMC) simulation. 
However, although certainly very stimulating, we find that these results are not conclusive be- 
cause of the following reasons: 

- PIMC simulations only give access to a fictitious classical dynamics of a system of polymers 
that represent the imaginary time evolution of the quantum particles. Despite some claims 
that this dynamics might be close to the true Schrodinger dynamics of quantum particles 
(see e.g. [22]), the connection between the two has never been convincingly shown. 

- We don't find the evidence reported in Ref. [6] strong enough to support the existence of a 
glass phase, even for the PIMC polymers. Our point is that, by analogy with simpler classical 
systems, in order to establish the presence of a glass phase one has to perform a careful 
systematic study of time-dependent density correlations, and show that such correlations 
develop a plateau at intermediate time scales, that there is a clear separation of time scales 
between a short time (so-called /3) and long time (so-called a) relaxation, and that there is 
aging in the glass phase. In absence of these evidences, avoided or interrupted nucleation of 
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As already stressed, our results for the glass (purple) transition line are reliable 
only in the high density and high temperature regime in the (p, T) phase diagram, 
where a semiclassical expansion makes sense. Since our theory is unable to in- 
clude exchange correctly, we cannot describe the interplay between superfluidity 
and glassiness taking place at low density and low temperature, say, within the 
dashed black line in the phase diagram of Fig. 2. In the rest of this section we 
describe two possible scenarii for the transitions taking place in the region corre- 
sponding to the question mark. 



3.1 A first order superfluid transition: the heterogeneous superglass phase 

In Fig. 3 (left panel) we report a conjectural phase diagram for metastable Helium 
4 in the {p,T) plane, as it emerges by the combination of Refs. [22, 23, 24] and 
the present study. The purple line is the glass transition line obtained in the present 
work within the first order semiclassical expansion. Below it, the system is a glass. 
In [23] it was predicted that the A -line might become a first order transition line 
at high density. Although this prediction has been made for a quite abstract lattice 
mean field model of dense liquids, it is interesting to discuss it since these models 
in several cases display a phenomenology that is very similar to particle systems. 
One should then observe a superfluid (SF)-glass (G) coexistence, leading to an 
heterogeneous superglass where regions of low-density superfluid would coexist 
with regions of high-density glass. The reader should keep in mind that the qual- 
itative behavior (in particular the form) of the A -line in the metastable phase is 
entirely conjectural. Moreover we currently don't have any estimate of the exact 
location of this line; therefore, its position in Fig. 3 is arbitrary. 



3.2 The homogeneous superglass phase 

A second scenario is that the A -line always continues to be a second order transi- 
tion in the metastable part of the phase diagram. In this case, it will very probably 
cross the glass transition line at some point, leading to a superglass phase (SG). 
This picture, which is displayed in Fig. 3 (right panel), is consistent with two quite 
different physical scenarios for the glass transition line at low temperature: 

- In Refs. [39, 40], a phase diagram like the one in Fig. 3 (right panel) was 
found. However, in the models investigated there, the glass transition is con- 
tinuous, like for spin glasses: the amplitude of frozen-in density fluctuations 
changes continuously at the transition. This is very different from the transi- 
tion of structural glasses, that we actually find in our semiclassical analysis, 
and that corresponds to a discontinuous jump of the amplitude of the frozen- 
in density fluctuations, the transition remaining however of second order. To 

a polycrystalline phase can easily be confused with the formation of a true glass, as it has 
been shown in a number of studies on classical systems. An example of this in the context 
of distinguishable Helium 4 particles was given in Ref. [15]. 

For this reason, the system of Ref. [6] might be in reality a moderately viscous supercooled 
liquid, and not a really arrested glass. We hope that further numerical work will clarify this 
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Fig. 3 (Color online) {Left panel) A conjectural phase diagram for metastable Helium 4 as it 
emerges from the combination of Refs. [22, 23, 24] and the present study. The stable part of the 
phase diagram is as in Fig. 2. The metastable phases are indicated in italic. The purple line is the 
glass transition line obtained in the present work within the first order semiclassical expansion. 
The 1-line becomes, at some point inside the supercooled liquid, a first order transition line [23]. 
One should therefore observe a superfluid (SF)-glass (G) coexistence. Note that we currently do 
not have any quantitative estimate of the behavior of the A -line inside the metastable phase, 
which is drawn at arbitrary position in the figure. (Right panel) Another possible phase diagram 
for metastable Helium 4. Same colors and remarks as in the caption of Fig. 3. The A -line always 
stays a second order transition. When it crosses the glass transition line, a superglass (SG) phase 
appears. 



match these two findings one can conjecture that below some temperature the 
glass transition line ceases to be discontinuous and becomes continuous in- 
stead. If this is the case, at low temperature glassy amorphous order would 
grow gradually from the liquid, and the relaxational dynamics approaching the 
glass transition would be very different from the one of classical supercooled 
liquids. 
- The second possibility is that the glass transition line continues to be discon- 
tinuous down to r = 0. This would imply the existence of a quantum discon- 
tinuous glass transition as a function of density. This has indeed been found 
in [21], based on a variational Jastrow description of the ground state of an 
interacting quantum particle system. Still some conceptual problems remain: 
indeed, all mean field models that could be solved exactly and that display 
such a transition in the classical limit, show that the latter becomes first or- 
der at r = [16, 17, 18, 19, 20]. The model of [21] displays such a transition, 
but its ground state wavefunction has only short-ranged two-body correlations, 
and therefore can be mapped onto a classical liquid; this is likely not to be the 
case for more generic models. A better understanding of perturbation theory 
around Jastrow wavefunctions, which was only attempted in [21], would be 
very helpful. 



4 Discussion 



Before concluding the paper, we critically examine the assumptions on which our 
method was based: we highlight its weaknesses and we discuss how to improve it. 
The first important issue that has to be discussed is the coupling between repli- 
cas. In section 2, Eq. (2) we have coupled the replicas as in Ref. [28], which 
is appropriate in the classical case. Here we discuss whether this way of cou- 
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pling replicas remains correct in the quantum case. In order to discuss this issue 
we go back to the original derivation of Monasson [26]. Since the glass is char- 
acterized by an inhomogeneous density profile, Monasson considered a system 
whose density is coupled to a static reference density field o{r). We denote by 
p{r,{x(T)}) = T.if{r — Xi{z)) a smoothed density profile at imaginary time T, 
where the function /(r) is normalized to 1 and very short ranged. Then in the path 
integral representation we get 



Z[o] = /"^{jc}e-^o[Wl-^^^''''^[''W-P('''W^)»l' 



(24) 



where 5'o[{ji:}] is the unperturbed action and the trajectories x,- are assumed to be 
periodic in t — ;> T + j3; a summation over the permutations due to Bose statistics 
is implicit. Following [26], we consider the partition function of the field a with 
"Hamiltonian" F[a] = —TlogZ[a] at inverse temperature j3m, given by Z„, = 

/ ^a{r) exp [ - j3mF[a]] = / S>a{r)Z[a]'". Integrating over a{r) we get: 

^ra.-I.aSol{^}]-j^„fdTdT'dr^,,[p{n{^{m-p{n{J'{T')})]^ ^ ^25) 

Defining F{r — r') = Jdzf{r — z)f{z — r') (note that F{r) is still a normalized 
short-range function), the coupling term has the form 



dab{TX) = j dr[p{r-{x''{T)})-p{r-{x\T')})f 

= £[FK(T)-X«(T))+F(xf(T')-x5(T'))-2FK(T)-4(T'))] , ^^^^ 



This coupling makes the configuration ^^(t) similar to JCy (t') or to a permutation 
of it. Thus, in the strong coupling limit, the particles build "molecules" made of 
one atom of each replica. If we make the approximation of ignoring permutations 
of the atoms of one replica belonging to different molecules, or in other word if we 
only allow permutation of entire molecules (molecules are Bosons but individual 
atoms cannot be exchanged), then we can re-label the atoms in such a way that 
molecule I is built by jcf , and so on. In (26), due to the short range of F{r) we can 
neglect all terms that couple atoms with ;' ^ j and we get 

datixX) ^ 2'£[F(0)-F(x1{t) -xf(T'))] ^ -£f"(0)(4(t) -xf (t'))' , 

where we used that F(r) is short-ranged and the coupling is large, so that atoms in 
a molecule will vibrate around the bottom of the potential well defined by F{r), 
and we can assume that xf (t) — xf (t') is small; the result is that the coupling in 
a molecule can be assumed to be harmonic (note that F"{0) < 0). Inserting this 
result in the path integral (25), we almost recover the path integral of the molec- 
ular system described by the Hamiltonian (2), except for the fact that replicas are 
coupled at different times. We therefore have to make the further approximation 
of ignoring the coupling of replicas at different imaginary times. 

The above derivation of our starting point, Eq. (2), is then based on two crucial 
approximations: 
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- We neglected exchange of atoms between different molecules. This is incor- 
rect when exchange is important. But taking into account these processes, the 
molecules loose their identity. A coupling term such as the one in (2) cannot be 
written since it explicitly breaks permutation symmetry of each replica. There- 
fore one should write a coupling that is manifestly invariant under permutation 
of the atoms of each replica, independently of the other replicas. This can be 
done, but solving the Hamiltonian looks a much more difficult task. 

- We neglected the non-local coupling of replicas in imaginary time. It is pos- 
sible than since the limit £ — ;> must be taken in the end, the precise form of 
the replica coupling is not relevant. It would be nice to perform the calculation 
including the non-local coupling, but this is technically more involved since a 
Hamiltonian formulation is not possible in this case. One has to work directly 
with path integrals. We attempted to perform such a computation but did not 
succeed. 

We note that although these problems seem very serious, they should not matter 
in the semiclassical regime where exchange can be neglected, and quantum fluc- 
tuations are small, so that the path integral is over a short time and coupling the 
replicas at the same time or at different times should not change the result a lot. 
Therefore, we hope that Eq. (2) is justified as the starting point of our semiclassical 
calculation. 

Another major source of concern is that a kinetic energy contribution to the 
equilibrium complexity appears in Eq. (15). This contribution obviously vanishes 
in the classical case, and luckily enough it still vanishes at first order in the semi- 
classical expansion. However, it is definitely non-zero beyond this order. This is 
problematic since it would imply that the kinetic energy of the liquid and the 
glassy states visited at m = 1 are not the same, which is instead a basic premise 
of the replica approach to glasses (all correlation functions should actually be 
the same). The issue is not the replica approach since the original Monasson's 
derivation can be formally generalized to the quantum case by using the density 
field p{r; {x{z)}) without encountering any difficulty, as we discussed above. The 
problem instead relies in the hypothesis behind the cage expansion, in particular 
the harmonic approximation. It is not clear to us how to find a solution. Possibly, 
a more refined description of the liquid and a more accurate expansion in e should 
at least reduce the problem quantitatively. 



5 Conclusions 

We presented a first attempt to build a replica theory of quantum glasses gener- 
alizing the recipe of [28] to the quantum case, and performing a semiclassical 
expansion that seems to give meaningful results at least for high density. In par- 
ticular, we estimated the glass transition line and found that, as already reported 
in [22, 23, 24, 25], quantum fluctuations promote glass formation. We critically 
discussed some issues related to our approach, in particular the role of exchange 
and a technical problem related to the kinetic energy. Although a full theory of 
quantum glasses must include these effects to provide meaningful results, they are 
not expected to be important in the semiclassical regime we focused on. 
Based on the combination of the present study and previous ones [22, 23], we 
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conjectured a phase diagram of dense metastable Helium 4 which is reported in 
Fig. 3 (left panel). In this case, it might well be that the glass ceases to exist in the 
strongly quantum region, since it undergoes a first order transition to the super- 
fluid; then, the present treatment could be appropriate to describe the glass, and 
the problem would be then to estimate the free energy of the dense metastable liq- 
uid in order to look for a first order transition. We discussed in Fig. 3 (right panel) 
another possible scenario, in which the glass phase becomes a superglass and a 
zero-temperature glass transition exist. The nature of this quantum glass transition 
remains unclear. 

There are several ways to go beyond our treatment. 

- One could try to use influence functional methods, which were used in [22] to 
compute the thermodynamical properties of the liquid. Although these meth- 
ods do not include exchange, they should allow for a more accurate description 
of the glass. 

- Another possibility would be to use finite-temperature variational methods [41] 
in order to extend the computation of [21] to finite temperature. This allows 
one to include exchange, but not to describe the normal phase correctly. 

- A careful investigation of the perturbation theory around the special model 
of [21] should help to understand if a quantum discontinuous glass transition 
indeed takes place. 

However, since the physical picture is very far from being settled, before pushing 
forward the theory it would be extremely interesting to perform Quantum Monte 
Carlo simulations or experiments (along the lines of [38]) in the dense metastable 
liquid region, to investigate the behavior of the A -line at high density. 

Acknowledgements We warmly thank S.Balibar, B.Clark, L.Foini, D.Reichman, G.Semerjian 
and M.Tarzia for many useful discussions. 
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